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ABS PRACH 


The topic of this paper is the stress analysis of an 
icebreaker bow. The purpose in writing this paper was to 
familiarize myself with the various methods that have been 
aeveloped in determining bow loads in ice breaking and the 
methods that have been developed to compute the stresses in 
a ship's structure. From a study of many methods, I wanted 
to select one that would be most useful to me in my day to 
day work when I depart from M.I.T. I selected as tne ship 
to model the United States Coast Guard Icebreaker WESTWIND. 

itiethie Pest SeClLlon. Ehe Vanrtous: metiods sot 
determining bow loads in ice breaking was studied. The 
various methods that were put forth by scholars and engineers 
over the past half-century gave bow loads ranging from one 
hundred and fifty-four tons to values as high as twenty-nine 
hundred and fifty tons of force at the pow in a vertical 
direction. I chose as the best results to use in loading 
the bow, those values which are determinable from the 
computer program developed by Dr. White in the doctoral 
thesis, "The Dynamically Developed Force at the Bow of an 

> Icebreaker." The computer program developned by Dr. White 
gives the vertical and horizontal forces acting on the bow 
as function of the length between perpendiculars, the beam 
at the waterline, mean draft, displacement, bow angle, 
spread angle complement, impact velocity, water plane 
@eemercient, longitudinal position of the center of floata- 
mrommetne Longitudinal position of the center of gravity, the 
height of the center of gravity, height of the thrust line 
above the baseline, bollard thrust that would be obtained 
for RPM being used during the crushing phase of the ice, 
lemgrtudinal metercentric height, kenetic friction coefficient 
~of ice and the compressive failure stress of ice. The verti- 
cal load determined by his computer program for the WESTWIND 
ramming solid ice at fifteen knots gave a peak value of two 
thousand tons. 


IPS ARCHIVE 
BEd 
peeved 2s e 






LIGRakY 
ta... POSTGRADUATE SCHOOL, 
rey, CALIF, 93940 DUDLEY KNOX LIBRARY 


NAVAL POSTGRADUATE 
SCH 
MONTEREY, CA 93943-5104 = 


The next section of this paper covers the various 
methods available for determining the stress loading inside 
the structure. The most used computer programs that were 
reported upon were the STRESS and FRAN programs. BoOcn OL 
these programs being strictly a frame analysis of the 
structure leave much to be desired in the complete know- 
ledge of the stresses present in a Complicated as UCtuLe 
Silene as an icebreaker. 

The most promising and most exact analysis that 
is available at the present time, utilizes the finite 

~element method. A program developed and presently being 
revised by NASA that quickly and efficiently analyzes almost 
mae structure is the "Structural Analysis and Matrix Inter- 
pretive System" (SAMIS). The remainder of the paper gives 
a short summary of the program, the data required by the 
emeaivyst, and finally the necessary steps required for the 
analysis of the icebreaker WESTWIND. 
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ii RODUCTION 


in the Ssareberider by von Ruskin, 89-1900, he 
Seates, "Take it all in all a ship of the line is a most 
honorable thing that man as a gregarious animal has ever 
Proauced. into that he has put as much of his human patience, 
Semmonm sense, forethought, experimental philosophy, self 
control, habits of order and obedience, thoroughly wrought 
Mandiwork, definance of brute elements, careless courage, 
careful patriotism, and calm expectation of the judgement of 
God as can be put into a space 300 feet long and forty feet 
Iigoaa," (Ref. 1)* 

Since that time, many ships have made history-- 
some of them by going to the bottom of the sea; others, in 
the defiance of the brute elements have eee eanesen the 
Arctic waste. This type of ship is what this thesis is all 
aout . 

The icebreaker hull at its birth was merely a 
standard hull with extra reinforcement in the bow. As time 
progressed, changes in design were based largely on the 
success and failure of past designs. Efficient ice breaking 
(or what was believed to be efficient ice breaking) led to 


the thirty degree angle bow. Removal of the unsuccessful 


a 
References are listed at the end of this report. 





bow propeller on the wind class breaker produced the step 
which prevented the vessel from riding completely up on the 
ice. With the advent of the computer, more efficient and 
better designs have been produced. 

Although very few papers have been written on the 
subject of icebreaker design, those that have been produced 
have made far reaching effects on the design of the ice- 


breaker hull. 
METHODS OF DETERMINING BOW LOADS 


The earlier works were done by the Russians. 
fee. karshis in his paper, "Ice Loads Acting on Ships,” 
(Ref. 2) put forth a formula to determine the impact load on 
the vessel. This ne considers the speed, ,the angle of 
blow, and the square root of the relative mass and the 
relative rigidity of the ship. He uses as an example a 655 
fom diSplacement vessel in contact with an ice floe that is 
25 meters in diameter and one meter in thickness. From his 
formula he deduces that the impact load is 220 tons. He 
assumed in this work a crushing strength of ice of 570 pounds 
per square inch. 

The general idea of this paper is that the speed 
of the vessel and the angle of blow in the area of contact 
with the ice floe are the major contributors to the impact 
load on the vessel; hence, if you double the speed of the 


vessel you double the impact load on the bow of the vessel. 


eine paper has a high theoretical content and requires many 
parameters difficult to determine. The relative angles of 
the hull form at the point of impact are necessary and these 
eam vary quite appreciably in a relatively short distance." 
feet. 2, pg. 5) 

ie NOG Lin cemapeis, | Impact Of  Satpe. vile 


fee, (Ref. 3) attempts to determine the reduction of speed 
o£ the icebreaker as it contacts the ice. An interesting 
sidelight of this paper is that the theoretical amount of 
force required to initiate the crack in the ice is quite 
pool, but that the force required EO _DEOpOgate tha crack 
might be quite large. It also goes on to show that the 
load required to break off pieces of the ice is more than 
three times the load required to initiate a crack. 

In his paper, Mr. Nogid gives a method to determine 
the speed reduction a vessel will have when coming into con- 
tact with an ice floe of given diameter if you know the 
maximum load the ship could withstand at the point of impact 
and the angle formed with the ice at the point of contact. 
The area can also be determined if you know what the strengtn 
of the ice is in the area in which the ice breaking operation 
is being carried out. 

He divides the forces and the strength requirements 
of the hull into two different areas. One 1S impact and the 


other is compression. 








Under the impact force, ne determines that this 
force is equal to the mass of the ship times the velocity of 
the ship, modified by parameters which are dependent upon the 
angle of the hull at the point of impact, the area of impact 
and other similar items. 

Under compression, the author assumes that the 
loading is directly related to the crushing strength of the 
ice. He also suggests some relationships between the changes 
of length and the strength of the side frames. His formula 
Suggests that the strength of side frames is related to the 
cube root of the ratio prototype to the proposed length. He 
also goes further in an attempt to relate hull strength to 
the frame spacing of the vessel. This paper, in general, is 
the theoretical approach to the old idea of taking a success- 
ful design and changing the dimensions to make a more capable 
design. 

U. N. Raskin in his paper, "Method of Determining 
the Stresses in Decks and Transverse Bulkheads Caused by Ice 
Loads," (Ref. 4) uses a compressive ice load on the side 
framing of a ship to determine what the stresses are in the 
decks and bulkheads. He divides the decks into strips with 
the only load carrying area being the deck itself and those 
beams in direct contact with the deck. The stringer is 


considered to be in simple compression, while the remaining 





wer ins act as rigid members to bending and carry only shear 
loads. In between these strips are imaginary elastic bands 
which allow relative movement back and forth between the 
bands. The bulkheads are considered as simple plates in 
end compression. 

Using these models, he presents a method to 
determine stresses in these decks and bulkheads. I believe 
it would be rather difficult to accurately model a complicated 
structure, such as an icebreaker in this form and achieve 
aeeurate results. 

In September of 1965, R. M. White completed his 
doctoral thesis, "Dynamically Developed Force at the Bow of 
an Icebreaker." (Ref. 5) In this thesis, Dr. White developed 
a computer program for the prediction of the dynamically 
developed force under the bow of an icebreaker while ramming 
Solid ice. 

The solution is based primarily on Newton's Law of 
Moecion. The SEepien was broken down into two basic phases. 
The crushing phase represents the local crushing of the ice 
to accommodate the bow. The sliding phase represents the 
Sromng Up of the bow without further crushing. "The final 
state represents (temporary) equilibrium when motion has 
See@pped; the vertical force at the bow at this state is 
relatively sustained and is the most effective in breaking 


fe ace. (Ref. 5 ig a3) 
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Dr. White gives one of the earlier developed force 


equations which is attributed to R. Ruineberg in 1888 (Ref. 6). 


= _ Tp (cos i, Cos B - - sin iy) 
BZ Sin 1, cos ocr fy cos 1, 
where 

Fay = Downward force 
T3 = Thrust available for ice breaking. 
Ey = Coefficient of friction (suggested as .05). 
B = The angle with respect to the center line plane 

of a normal to the shell. 
ly = The angle the stem makes with the base plane. 

Dr. Wnite states that Ruineberg assumed no momentum 


effects and no forward motion, so that all thrust is applied 
to ice breaking and is applied horizontally at all times. MThe 
direction of friction force remains the same during forward 
horizontal progress and the trim does not effect the solution. 
For the icebreaker WESTWIND, this would give a downward force 
eeeerour hundred and fifty-four tons for this continuous mode 
eeeuee breaking. 

The next equation Dr. White discussed was one 
developed by A. Kari in 1921 (Ref. 7). 


7 — 4480 ACL sin 8 
Bz H 





where 


A is displacement in tons 

O aserChandger ine ceim 

L is the length between perpendiculars 
H is the draft in feet 


and 


where 


BM, Soeassimed sto penwcaual Lo GM, . 
GM is longitudinal metacentric height in feet. 


Kari makes the following assumptions: (1) no 
momentum effect, (2) the bow rises the thickness of the ice, 
(Momene distance from point of contact to the center of 
Flotation and the center of gravity are the same, (4) 
effective displacement and the draft is not changed by load 
force, (5) the center of flotation and metacenter remain 
mie, (6) C = .07, and (7) there is no frictional force: 
This would give a downward force of approximately seven 
nundred and fifty-eight for the icebreaker WESTWIND. 

The third method discussed by Dr. White was one 


womerby F. R. Simonson in 1936. (Ref. &) 


= 
ibis: 
F om 


BZ tan (iB + 9) 





In this method, Simonson assumed there are no momentum effects, 
fection 1s negligible, thrust is directed horizontally, the 
coefficient of friction serves as a pivot point, and there is 
no change in displacement. Dr. White states, "... Simonson's 
equation is limited to being a good approximation for the 
Seopeed Equilibrium position." (Ref. 5, pg. 46) This would 
give a downward force of one hundred and seventy-three tons 
for the icebreaker WESTWIND. 

The fourth work discussed by Dr. White on the 
downward force developed in Lce breaking 1S one by Le V- 
Vinogradov which was published in 1946. (Ref. 9) Dr. White 
States that this is the first time that force due to ramming 


was put into useable mathematical form. 


** 
a ae ve [l - (1 - e* sin*i,) - v4] 
ee Sh te Lx eee , Gh 
where 
fe 
_ cea 8 tan 1p 
x. = 
“he 
Lee cot 1, 
ite tlic! ail ce O1ls 
Y=- - 
She 
1+ ee COL lp 


K* 
Symbols have been changed to previously defined symbols. 





W = displacement 


D = GQrett wn oft: 
(e k 
Bec lte we) 
W 2 W 
Cah = block coefficient 
Cy = waterline coefficient 
bs q 
ae) OD 
q -—mradtstance £rom polmte Of Jmpact Eo Center of flotation 
= —— cM -H- Cc 
foe. CL “L B 
W 
ee = speed prior to impact 
‘oe speed while sliding up 


Vinogradov assumed that thrust is always horizontal, 
change in trim and draft is small enough so as not to affect 
the waterplane characteristics and the metacentric height, 
and GM, = BM. 

In Vinogradov's equation the force downward is 
related to the frictional force, the angle of the bow, the 
spread angle of the bow, the coefficient of friction, the 
block coefficient, the waterline coefficient, the distance 
of impact from the center of flotation, the speed prior to 
impact, and the speed while sliding up. This would give a 
downward force of approximately two thousand nine hundred 


and sixty tons for the icebreaker WESTWIND. 





The fifth equation discussed by Dr. White is one 
Pyert. Richardson in 1959 in some personal correspondence to 
Dr. White. "The development was almost identical to 
Vinogradov's but did modify some of his weaknesses to some 
extent. For example, Richardson uses a term for the loss of 
energy due to wave and frictional resistance (not ice) from 
the instant of contact up to the moment the ice breaks or 
motion ceases. He also recognizes an effective increase in 
ilwmiass Of the 1cebreaker due to entrained water.” (Ref 3p 49) 
The final equation discussed by Dr. White was one 
developed by V. R. Milano in 1962 which is a modification of 
Vinogradov's equation. Dr. White states, “One of the main 
Pemtributions was to express thrust as a function of ‘Bollard 
Pe@eeoe°§6©6(ReE. U5,lhUp. 6S) 
Dr. White in his analysis assumes: 
1. The force normal to the plating is 
represented by the product of the area 
of contact and the compressive failure 
stress of the ice. 
2.  -wnere 1S a. -LPiCction £OrCe Acting ps nseme 
plane Of Ene plating. 
3. The icebreaker is treated as a "solid body.” 
4, The radius of gyration of an icebreaker 


Gan be-assidqned 0.260. 





5. Crushing has ceased when velocity is in 
the direction of the angle of the bow 
eR Se hice NG liomOt uti 11m, 

Gi. Te sliding phase, the point of 
contact is fixed and is at the level 
of the waterline. 

Results show that the peak load occurs during the crushing 
phase. 
In his computer program, the downward force is 


determined as a function of 


Computer 
Symbol 
L Length between perpendiculars, ft. BP 
B Beam at waterline, ft. B 
H Mean Giant, fT. H 
A Displacement, lbs. DIS 
lp Bow angle (from base line to stem), radians BA 
SA Spread angle complement (normal to bow SA 
plating with respect to center line 
plane), radians 
Va Iinmpace velocity, £t/sec. V1 
a= Cy Oo, Waterplane coefficient, dimensionless AL 
LCF Longitudinal position of the center Of Gr 


PlOtacton = sit “AbcoOr tamMigsniOS>. a ans 


EOEWara), ft. 


ee ial ts 





Computer 


Syimbow 
CG LCG, Longitudinal position of the center CG 
CieGpay ley (nk abc Ob amladchips ~ tie 
Ie pennrsiel@:)) 6 ikea 
KG Height of center of gravity above base GK 
Heine. Ee. 
d Height of thrust line above base line D 
near center of gravity, f£t. 
TROL, Bollard thrust which would be obtained TB 
for rpm used during crushing and 
Siicuna, Los, 
GM, Longitudinal metacentric height, ft. GM 
fy. Ronemve COCLEIC Tent Of EricblenwE FK 
1ee/ ship. 
Oo Compressive failure stress of ice, H/Et. SIG 
The following pages are Dr. White's complete 
Peegeam for developing this downward force. (Ref. 5, 
moemeso5-393). 
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PRE oo ANALYSIS OF SHIPS STRUCTURES 


In the analysis of the Coast Guard icebreaker 
SEAclER, in Lloyd's Register of Shipping Research and 
Technical Advisory Services Report No. 5095 (Ref.10), the 
IBM program STRESS was utilized. This is a computer program 
Pome performs a linear elastic analysis of a framed struc— 
ture under the influence of a static load and is capable of 
eieing two or three dimensional problems. 

The condition of loading that was utilized was 
400 psi ice pressure on the midship section of the icebreaker, 
and the ice condition was considered to be 10 feet in thick- 
ness. There were two conditions imposed with this 10 foot 
thickness of ice. 

(1) Two feet of ice extending above the waterline 
and eight feet below it. 

(2) Assume that all the ice was placed at the weakest 
pelnt of the structure that was likely to encounter 
ice load in the ice breaking operation. 

Once again in this program, you have a problem of 
the plating not carrying any load which certainly, in the 
analysis of an icebreaker where you have 2 1/2 and 3 inches 
@emerating, will lead to large errors. 

in the analysis that: was Carried owt, Cwo Sees oF 


end constraints were considered: 





(1) Supports only at the extreme corners of the 
SeruUctire Or Vor nave FOUr JOllltse. 

w2) Supports at all inboard ends ef kKeeteens, crc. 
ineadaLteion to vthe Condition abvovergoind coma 
Speci yOumks. 

The bending moments obtained from this program with 
four and eight supports varied markedly with factors as low 
as two and with others as hign as 60 and none of the varia- 
tions showed any real relationships. The four support condi- 
tions gave values at various positions in the structure that 
at some points were much greater and sometimes much less than 
mirewcOorcresponding points of the eight support conditions. 

The results of the paper did not give any real 
concrete evidence as to what the actual stresses would be in 
the vessel considered under the 400 psi ice load condition. 
The conclusions were only that the structure should be of 
the grid structure instead of the truss structure. Once 
again, this is nothing more than a selection from the best 
of the successful designs to be applied to a new design. 

Concurrent with the previous report, further 
investigation was also being done on these amidship frame 
areas by Paris Genalis. (Ref. isl, FFs also used the IBM 
computer program STRESS to analyze the amidship's section. 


Ree uOwes ein Adc ewe pO: eeem 


Im Mis report, he breaks the analysis down into 
mouUm Stages: 
a) Consideration ef, part Of only oneerrame. 
(2) Consuderatweon of one half of one frame- 
(3) Two dimensional analysis of one complete ring 

Frame. 

(4) Three dimensional analysis of more than one 
complete ring frame. 

Bach of the analyses carried out was used to 
substantiate the assumptions and methods utilized in the 
next analysis. In the final three dimensional analysis, the 
plating was taken into account by replacing it with stiffeners 
to get away from the problems that Lloyd's of London had in 
mierc report. 

Each of the problems considered in the three 
Gimensional analysis, took a total of one hour or more of 
computer time. The number of calculations was very large, 
hence the Selanne of error was very large. 

Mr. Genalis also went further in his analysis in 
that he considered machinery weight, steel weight, and the 
effects of varying amounts of buoyancy in addition to the 
ice loads. The wide variations in the results of his analy- 
sis point out the absolute need for proper modeling of the 


Maate rice in order to determine realistic Stresses. 





Further limitations of the STRESS seregram 1s Ehat 
it cannot handle curved members and thus assumes that each 
member is straight and slender. 

Another application of the IBM program STRESS was 
completed by Lloyd's Register of Shipping, Research and 
hmeemircal Advisory Service Department in Report No. 5051, 
Meet. 12), in which a structure of shell plating subjected 
to ice pressure was analyzed. "The purpose of this investi- 
gation was to design a simple panel of shell plating stiffened 
in such a way so that the overall and local strength was 
sufficient to withstand a specified ice pressure uniformly 
distributed over the whole pvanel and for a concentrated load 
acting on one frame only for the full depth of the panel." 
Ger. 12, p. 1} 

Here again, we have results that do not take into 
account the shell plating other than to simulate it as part 
of the stiffeners. It does give results that comvare well 
with simple pean theory, but does not give the information 
necessary for the designers who must, in these days of 
economy, depart from the large factor of safety of five or 
ten to insure that the vessel will be successful and safe 
from the environment it operates in. 

The IBM STRESS program was also utilized by 
@entsaiter, inc. of 1725 K Street, NsW., Washington, Daye: 
in "Polar Icebreaker Preliminary Structural Design and 


Seecial Studies,” (Ref. 13) completed in August of 1968, 
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and gave a summary of known ice properties and an analysis 
Of possible structures to be utilized in a new type of 
icebreaker called the M-10 at that time. 

The ice loading used was that developed from the 
pimeeivysis by M. K. Tarshis, “Ice Loads Acting on Ships, ~ 
free. 14) a translation of a Russian text published in 
Powe cransport, Vol. 16, No. 12, 1957, pg. 19. 

Using these values of ice loads they determined the 
capabilities of the proposed design along both elastic limits 
maemolastic hinge limits. They included factors of 1 1/2 £or 
modest loading and two for static condition loading to 
Beecunt for the lack of reliability of the values of the 
properties of ice. 

The results once again lack aspects of reality due 
Gemene limitations of the STRESS program. 

During recent years a new method called the finite 
element technique has been developed to give more accurate 
solutions for structure analysis. The first easily under- 
stood, comprehensive, presentation of the method was written 
by 0. C. Zienkiewicz in his boom, "The Finite #Hlement Method 
in Structural and Continuum Mechanics." (Ref. 15) Develop- 
ment of a computer program capable of analyzing a wide 
variety of structures was developed by the Jet Propulsion 
Laboratory, California Institute of Technology, Pasadena, 
@alifornia, called the Structural Analysis and Matrix 


Interpretive System (SAMIS). (Ref. 16) 
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THE STRUCTURAL ANALYSIS AND MATRIX INTERPRETIVE SYSTEM 


The Structural Analysis and Matrix Interpretive 
System (Ref. 16) uses elements that are at present restricted 
to flat triangular facets that are joined together along 
their edges,and line elements that are joined to the rest of 
the structure at their ends. All elements are capable of 
resisting stretching, shearing, bending and twisting loads. 
Heating, acceleration and pressure loads can also be analyzed. 
Additional loads can be introduced as energy equivalent con- 
centrated loads at points on the structure. 

Structural changes are defined by stresses, deflec- 
tions, flexibilities and stiffnesses. It also has the capa- 
bilities of computing the natural frequencies of the struc- 
ture and the mode shapes. 


PHiertinite element anvolved Utila Zese eyo eta o. 


1. The structure is divided into small elements 
2. The problem is solved by a structural stiffness 
analysis. 


Load deflection relations are taken for each element of the 
structure. The coefficients of these relations form the 
stiffness matrix. When the whole system is taken into account, 
the load deformation relations for the entire system stiff- 
ness matrix is developed by summing the stiffness matrices 
SErthe pieces composing the system. Where there are common 


grid points, the forces are simply added. Boundary conditions 





are also formulated in matrix notation. The displacements 
at each of the corners of the facets and the ends of the 
beam that make up the complete structure are determined by 
solving systems of equations simultaneously. By taking 
these displacements, the stresses in each of the elements 
are determined. 

In the Jet Propulsion Laboratory Technical 
Memorandum No. 33-317, (Ref. 17), the mathematical equations 


of the structure are given as: 


[K] {a} - EST as Ci ic SSI ce) [M1 {d}{o} [Rl] {d} 


l| 


where {d} 


II 


The vector of grid point displacements 


(dOotS Indicate time derivatives). 


Gi = The vector of element stresses. 

[K] = The small deflection stiffness matrix. 
[C] = The small deflection damping matrix. 
[K.] = The initial stress stiffness matrix 

r = A scalar defining the magnitude of the 


IMLeVvalL! Stress Grctre wove ron. 


oO 
— 
= 

6) 

— 
I 


The mass loading matrix. 
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The-fowce loading Matix. 


Pe =/2Pne Marri Ol Stress CoOcrtrlcremes- 


Eoundary Conditions are applied on "d" and "A" to make the 


above equations solvable. 
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The static and dynamic displacement response of 
the structure is given by the first equation and the stress 
associated with a given set of displacements is given by the 
second equation. 

An assumed displacement function which is continuous 
along a facet's edge and continuous over the facet is used to 
@ewerop the stiffness coefficients by minimizing the potential 
energy of the entire structure. The resulting stiffness 
coefficients define the forces and moments at the apexes of 
the triangle which satisfy macroscopic force and moment 
Seu torlum conditions. 

The displacements used to develop the stiffness 
coefficients are also used to develop the loading coefficients 
under the restraints of minimum potential energy. 

The stress coefficients are developed using the 
same assumed displacements functions and stress-strain and 
strain deformation relationships. The stress so developed is 
the mean stress ate tNestacet Cenerord. 

Two coordinate systems are used to implement this 
program. The local system describes the grid points of each 
element as to their local relationship, and the common system 
describes the overall system. In the facet element, the local 
x-y plane is assumed to be coincident with the midplane of the 
facet. The displacements of each grid point are defined as a 
Mee vector which consists of translations in the x, vy, and Z 
epeeoctions, and the rotation about tthe x and y axis. Each 


component of these deflections is assumed to vary linearly 
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over the midplane surface. The deflection component at any 


point is therefore expressed as 


6) ees oe a ae 81 + Cc. 
aL a Ly rf 


iecie stress-strain relation, it is assumed that the stress 
in the z direction 1s zero; assuming that the thickness 
radius ratio is small compared to one for the shell. 

The selection of the triangular elements must be 
such that the stiffness matrix is positive definite. This can 
be done as long as the largest angle for any triangle is less 
than 90 degrees. If an angle is larger than 90 degrees in the 
meet, the stiffness matrix will be indefinite and thus cannot 
be used. The most accurate results are obtained if the 
triangles are all equilateral. 

Discrete loads acting at a joint are described by 
forces acting in the x, y, and z directions and moments acting 
about the x and y axis. When torque is considered, it must 
be treated as a couple with forces in the x and y directions. 

Local coordinate axis are used in developing the 
Matrices of coefficient and, hence, the need for a trans- 
formation from the common axis. This 1S automatically 
accomplished if the grid points are given in the data as 
the common coordinates. 

The line elements used are superposition of models 


imemeaxial elongation, torsional rotation, shearing and bending 


Where stiffeners provide the resistance to bending, the 
classical bending element is used in the analysis. When 
using a facet, the shear bending element is used to provide 
the necessary stiffness from the supporting frames as this 
pEovides deformations which are consistent with those of 
Eme facet. 

As in the facet element, assumed displacements 
are used in the development of the stiffness and loading 
BPeemricients. The stress coefficients are obtained directly 
from the stiffness matrix. 

In the line element, the local x axis is considered 
to go down the centroid of the beam. The program automati- 
cally makes this transformation if the line element grid 
points are given in the common coordinates. 

Due to the versatility of the program, the input 
necessary from the analysis for the SAMIS porgram is 
e<tensive. The first section of input for the program is 
Ee program control cards called "pusedo instructions. ~ 
These cards control tape assignments, matrice maming and 
matrice manipulation. Card format is divided into alpha- 
numeric names of the matrices, tape assignment instructions, 
instruction sequence number and program control information. 

The card format is broken up nto ten f1elde. 
mae 'O" field is the instruction sequence number. The “Ay, 


"Bp" and "c" fields contain the tape storage assignments, 
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the fields numbered 1, 2, and 3 contain matrix names, the 
"code" field contains the subroutine or subprogram name and 
mice t field contains the control information for the 


operations to be performed. Card format is as follows: 


ieee. 15) 2X,A3,,103)3X%,15 | 2X,A3,3) 4X, A4 | 3X 2k, AS 3 ea 





ae Ey, RSS a ome nn. | 





Pemeqao instructions are: 


ee ae 





TABLE L 
(Ret 16,79. 144) 


OPERATION PSEupO INSTRUCTION. 


imastructions Interpretation 
ADDS EP OletieGo <6 af BZ 
BILD Construct small deflection stiffness, stress, 


loading, and/or mass matrices as Al, Bl, Cl, 
and Fl, respectively. 


CHIN Form B2 such that B2 Boe = Al and C3 = a2 
where Al is symmetric and positive definite. 

CHOL Form Gs = re B2 using Choleski decomposition. 

CONT Continuation Card (see pg 43) 

CODE Transform Al to coded format as C3 

COLS Put .the Al mater in column lvetrng andmea rm 
Hie Gos. 

DECO Transform Al to precoded format as C3. 

FILL Read Al, B2, and/or C3 anborcorce. 

Pie Form C3 = aie 

INKS Print Mathices only 2s and orre. 

MULT Form @€3. = Aa eZ 

READ Read matrices Al, B2, and/or €3 fron carcace 

ROWS Put the Al matrix in row lasting and call ae 

ROOT Find latent roots and vectors of Al, a 


SYMMECLIiC MatEix. Let vectors be B2, GOoOES tae 


SAVE WeitevA lye B2, and/or rCs on eape- 

SORT Sort a matrix Al a COW -OENCOLUMI Nasser 
SUBS POU CG k=) Atl eaeZ 

WASH WASH Al elements from B2 to produce C3. 





Psuedo instructions called logic instructions 
Peevide the Capability o£ loop operations and are also 
mequired for transfer of control to the success exit. When 
logic instructions are used, the "E" field indicates the 
number of times the logic instruction is to be carried out. 


iegie instructions are: 





TABLE 2 
(Ret 2 Bey ea 


BpOGIC PSEUDO INSTRUCTIONS 


iseructions Interpretation 
PREP Prepare for multiple execution of the 
following instructions. Execute the 


instructions between the PREP and the 
next, BACK instruction the’ nunmoerso. 
times specified in the "E" field. 


VARY Vary matrix or tape numbers in the next 
instruction by augmenting corresponding 
Field data by the specified integers 
after one pass. 


BACK Back up and repeat instructions after 
PRE 

RRS Disrupt errors can be corrected as 
follows. 

She Skip the next "n" pseudo instructions 
where "n" is specified in the "E" field. 


(Skip cannot be included between a PREP 
and DAGK INS truce uenoe. 


STOP StOp This Case and go) Lom EL. 
PAWS Pause in the calculations. Operator can 
restart at any time. (Not overative on 


TBM 7094-7040 DCS). 


HALT Halt and indicate “a successtulme re, 





Control cards needed to direct computations for the most 
general case of static and pseudostatic loading and normal 


modes 1S a single analysis would be: 





SAC HE 
9004 
9003 


200 1 


JIAEOR OME 


SAOOAE 


2G 


MEO'OZ 
HOO 1 
POO 2 
E2003 
| ARO ONS: 
11004 
AILS NO) 
eZ O0)2 


9001 


WERL 


KER2 


KERIL 


Ker 


DURL 


Reel go 


WTRL 


WURL 


KURI 


S100) 4 


10001 


OZ 


11001 


OZ 


12003 


9002 


11004 


F002 


12003 


OOS 


ARCI 


KURL 


WURL 


WLCL 


Pee 


VDC lL 


Veer 


DWC1 


AUC2 


FORTRAN STATEMENT 


Cor 


READ 


BILD 


ADDS 


500e 


12001 


12 Oe 


One) 


OG 


KERL 


Rai 


ATCL 


RBCIOL 


WLCL 
Tee 
Dect 


DDC1 


VCea 
DWC1L 
AUC2 


Drew 


-yo5 
le 


510 Ny 


SAG, 


7010 


890 


1 - Facsimile of Operational Control Cards 
(including normal modes) 


(Reps. 19% ODS ancy) 


L 


Pag. 


ae ae 





NOTE: 


HOO: 
200.2 
E2002 


2001 


70.0 1 


ZOOS 


11006 
12002 
2002 
LOOOL 
EO O19 


2 OO 


EOI 


EOC! 


KURL 


Bpew 


DDC1L 


FOC2 


BORZ 


WTRI 


Foes 


FOR3 


Cee 


DDCL 


DTe1 


SCR2 


Ree 2 


RDEZ 


FORTRAN STATEMENT 


9002 


12000) 


T2002 


2 OO 


2) 0) 7 


DECI 


BOCe 


PHeL 


BiLew 


Dac i 


Bret 


READ 


ADDS 


CHOL 


CONT 


COS 


Pie 


MULT 


MULT 


else 


MULT 


DECO 


INKS 


FILL 


MULT 


psc © 


INKS 


HALT 


coon 


9002 


IEEO OG 


LOR, 


eZ O03 


LUGO 


ZOO 


te OD 


12001 


FOC2 


DEC} 


GCC1 


FOR2 


GCC2 


INGE S 


FOR3 


GCC3 


DDC1L 


REEZ 


RDEZ2 


10 


305 


yz, 


eae 


Replace N by the number of structural elements and P 


loan NG ie Cran 


Replace N, by the number of cards of codes associated 
with prescribed displacements (12). 


eee 





Statement number 9, for example, means take 
fae tx WUR] om tape 11, position two and multiply it by 
matrix TEC1 on tape 9, position two and place the resulting 
Meet x On tape 12, position 3 and call it DCCl. Instruc= 
Exons are performed in the sequence that they are placed in 
the input. Any Similar type program can use the same 
sequence. Comment cards may be included in the pseudo 
misc ructions. 

The next section of input defines the model weight 
at each grid point which is desired, in addition to that 
automatically modeled by the SAMIS program from given data. 
The first digits of the row and column numbers correspond to 
the grid point number,and the last digit gives the corres- 
ponding degree of freedom the weight is associated with, i.e. 
ime eand 3 the x, y, and zg direction respectively; 4, 5, and 
6 the off center weights about the x, y, and z axis. This 
Mereiax 3S Used to introduce concentrated loadings. Card 
format is as olen 
Soca one 


on [w [eels [=] «[« [s 


Nom OF Matrix | Row(-1) on Precoded (1) or 
GSacas Orwlasize column(l1) jcoded (0) 
eu listed 











fara three 






ite BAG 9 


——_ i= oe 


fewrre values row or column Listed 


For coded format 


Card two 


Row Or column or Value or 
column and row rem eng il>.< 
listed element 

The next matrix data input are the nonzero 
accelerations associated with the restraints that are not 
modeled by the SAMIS program. Corresponding to the assigned 
loading conditions, gravity (g) accelerations are supplied 
[iene xX, y, and z direction, and (g)/in accelerations 
forces are supplied about the x, y, and z axis. This matrix 
is utilized to make equilibrium checks. Card format is the 
same as the weight matrix. 

The next section of input gives the nonzero 
elements of the stiffness matrix that are not part of the 
elements listed. This matrix gives the holonomic boundary 
conditions imposed. Card format is the same as for the 
weight matrix. 

Two cards giving the modulus of elasticity and the 
modulus of rigidity and a card of zeroes complete this 


section of input. 





fapulated grid point coordinates arervElcnegiven 
in ascending order. If the coordinates are given in the 
common coordinate system, the SAMIS program makes the 
necessary transformations to local coordinates for genera- 
wom OL the various matrices. One grid point is given on 


each card with the following format: 


i 
eed. NO. Cri DOM: Wer. 


Following this section of input is the element 
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Grid, poInt Coordinates 


data cards. This input describes each element as to its 
type, position, size, geometry and the material type it is. 


The card format for each of the elements is as follows: 
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€ € € 


ARTRUSPI Z : A De 
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*pr0o0p qUTOd PTAD *sqng paTuL *pz005 
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aoezans Seen 
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FACET ELEMENT DATA 


Format 


A Gk 
IES: 


ry 
~] 
© 


er aera 8, 


x 


* 
Leen 


it 


Sal 


50 


FIRST -GRzD 
POT 


TAT RD Gi Poi 


Interpretation 


Indicates, Gard Munser =a 
Element number. 0 < K o 999. 


Specifies the Facet element 
assumptions. The first digs, 
3, indicates three grid points 
are required. 


First ~ SeCond , alc send ac 
elastic grid soink numbercmeon 
the element. If negative, 
solution displacements are 
produced in the local cgere mie 
ate system at the Grid peri. 


Data di tis fle la ieee 
ignored. 


First, second, and third sub- 
stitute grid point numbers. 

Tf the grid point numbers are 
negative, solution displace- 
ments are produced in the local 
coordinate system at the grid 
POLE. 


Facet mass 5 
M 0, total mass:#sec /in 
M 0, mass per unit area: 
tsec2/in 


[VIA 


Bacet Enieckness in anenece 


Name of the structural maceria- 
The first Ewo characters on seme 
material name must match tne 
£LiyvySt EWO Characteramor gem 
material name in the material 
table. 


* e a e 
Momoers gavyen in this column are to be taken literally ae 
user must substitute appropriate numbers for letters given. 





TABLE > (Cone "d) 


PaACET ELEMENT DATA 


* 
Cara Columns Format item interpretation 
st sk 2 Indicates card number = 2. 
2-4 ita K Plement nomber. 0. <a 
(may be omitted) 
2 5-7 ee ie Sul Same as Card 1, columns 5-7. 
(may be omitted) 
eZ 8-14 tyra CA Continuity. boundan, come. 
tions at the first, secemna, 
foo 2s £720 CB and third gnpid  ooincs. Welt soem 
substitute points are given, 
Ee) 0 EC these apply to the substitute 
points. If only elastic#appea, 
they apply to the elastic. 
DS) as 0 == Ignored 
Bo-22 Ene P Normal pressure: pounds/inch” 


Positive in the plus z 
direct lon: 


2 43-49 Hi 20 si Temperature (Degrees Rankine) 
of the material (used to de- 
fine elast icuconsctantes 


2 50-56 ara, Ab Upper surface temperature 
change (Degrees Rankine) from 
Zen Stress —eCMpebaeurise: 


2 a7 65 area, T Lower surface temperature 
change (Degrees Rankine) from 
ZGeGO0 SELCSS eemperacicime 


2 64~-70 B00 ve Ignored 
2 eee 2 A2 ——- Ignored 
3 Omit 
4 Ore 
x 
Mumbers given in this column are to be taken literally. ine 


user must substitute appropriate numbers for letters given. 











EACH LEVEN) Dire 


* 
Gard Columns Format Item Interpretation 
1 Tl 5 Indicates card number = 5. 
2-4 16s: K Etenent numbers (055 ek zoe 
(may be omitted) 
5 5-7 as Sa. Same as Card 1, coiumns 5-7. 
(may be omitted) 
5 8—70 OG 7 20 Coordinates (inches) in the 


overall system of the origin of 
the local coordinate system 
(X1-Y¥1.21)? a point on the 
local x axis (X9,Y9,22); anda 
POlnE in the x-y piane (37 ae 
23) MOnNCOllinéGar svikenie de. 
first two points are selected 
to define desired direction of 
the Local’ za .axis. 


5 PA 72 A2 a Ignored 

6 Oni te 

ui di Loa 7 Indicates card number = 7. 

4 2-4 ib K Bbemen t Gitlin mre 0 < alae 999. 
(may be omitted) 

| 5-7 i Grae Ses Suk Same as Card i, columns s.7- 
(may be omitted) 

Zi 8-19 97.0 Coordinates” (2necics) tomeer 


substitute grid points corres— 
ponding to the first, vsecene™ 
and third elastic grid =peimess 


7 71-72 A2 G Coordinate identification. Le 
C =} ih, ®substitute guid pels 
coordinates are in the local 
System, Tt © AIL, coordi. 
ates are overall. 
8 Omit 
* 
Mmmbers Given in this column are to be taken literally a ime 


user must substitute appropriate numbers for letters given. 





* 
Card Columns Format Item Interpretation 
Tl S Card number = 9. 
= He. K Elemente number Uh ss yor 
= x Le el same as Cand i, columns S275 
(may be omitted) 
9 o7.0) 947.0 Coordinates (inches) for the 
First (xX7,V1,2]), Seeene (X90, 
V5) 29) ane Enc a op 
elastic grid points of the 
Facet in the overall or local 
System. 
g ei A2 E Coordinate identification.) se 
C = L elastic grid point ™ceor- 
dinates are in the local system. 
Tf C 7 i, -coerdinates oreo ea. 
elle. 
“Numbers given) in this column are to be taken Tvterat aie 


TAeEE >  VeCont a) 


PACE ELEMENT eibpe 


user must substitute appropriate numbers for letters given. 





Aer 


(Ret. L667 pas. 1235. ci 2s Oe leone) 


= LINE ELEMENT DATA 


ileal 
GAGE D, 


POINT Pd 
Ss 





SECOND GRID 


ea POINT 


ne 
Card Columns Format Ltem* Interpretation 
i 1 ak al Indicates card number = l. 
il 2-4 ES K Element number. 107 <i <eee 
il ja 1; eZ 2A Specifies the lineelement equa- 
tions. Rod and tuke equations 
and (A = 1) elementary beam or 
(A = 2) shear beam displace- 
ment assumptions. 
al: 8-14 7.0 GA First, second, gand: tiie gia 
point numbers for the element. 
Is) 211 E70 GB The grid number of the third 
qrid point may se omitreaei. 
Paes Fea GC numbers are negative, solution 


displacements are in the local 
coordinate system at the grid 
porate. 


ik 29-35 1 ale A Cross-sectional area (axial) 
in the y-z plane: 1neie=—- 


1 36-42 Ey 20 J Torsional rigidity against 
Ewast oO Olt tem ax ioe 
inches =<. 


1 43-49 ice a8) A Effective area deforming in 
shear in the x-v plane due to 
avy Loree ™anehese. 


il 50-36 7.0 iL Moment of inertia resisting a 
Moment abeubky tne) Z axis. 
THenes == 


ss : ’ 
Mumoers given in this column are to be taken literally eine 


Meer must Substitute appropriate numbers for yleerernc: 
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TABLE 6 (COrme’d) 


LINE ELEMENT DATA 


Gard Columns Format cen Interpretation 


ih 57-63 i 7.0 A BELECELVE Sheareetea = foie 
ap s s 
shear force tre etesZ dicee— 
tion: inches2, 


1 64-70 B30 a Moment of inertia resisting a 
¥ twist about: tie vy anas: 
inches4., 
Ih 71-72 A N Name o£ the structural mater- 


ial, The first two Characters 
of the name must match the 
first two characters of “ene 
material table material name. 


2 i ae Ae 2 Card number = 2. 
2-4 nS K Element numbers. 0 sk ee 
Z 5-7 XZ 2A Same as Card 1, columns 5-7. 
(may be omitted) 
2 B21) Zig 0 CA Continuity boundary condi tiems 


CB at the first and second grid 
points. LiMstesti tutes gare 
voints are given, these apply 
to the substitute pointe ae 
only elastic appear, they 
apply to the elastic. 


222.0 E7.0 = Ignored 


29-35 EO M Line element mass. If M < Q, is 
the total mass (#sec¢in). If 
M> 0, M is the mass per unit 
length (#sec*/in@). 


2 36-42 E7 .0 P Normal pressure: pounds/inch 
(positive in the plus z 
direction). 

7 SS, E7.0 Ty Temperature (degrees Rankine) 


o£ the material (used tovde— 
fine clastre. cons tanca = 


2 D020 ow eae) aE Mean temperature change (de- 
grees Rankine per unit of 
€ross-Sectilonal arcay— ye mo 
the element from tne zero 
stress temperature. 


* s 
Pmmoers Given in this column-are to be taken I[lterall, ane 
user must substitute appropriate numbers for letters. 
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GABEE 6. (Conta) 


LINE ELEMENT DATA 


Card Columns Format Ttem* Meer pee ra eon 


2 Byblos: yee 0) dl Temperature gradient (degrees 
¥ Rankine per Unlieeucss— 
sectional moment of inertia 
Ie! in sehe 7 diineec elon. 


2 64-70 B7.0 sh Temperature gradient (degrees 
Rankine per unit cross- 
sectional moment of inertia 
I) im the y diveeeion. 


2 Tilers Jig A. —— Ignored 

3 Ona t 

4 Omit 

5) i al: 5 Indicates card number = 5. 

5 2-4 13 K Element NumMpeme Oc ae 
(may be omitted) | 

5 5-7 a2 2A Same as mlards, Columnss 7s 
(may be omitted) 

S Sete, SHE Fe, == Coordinates (inches) in the 


overall system of the origin 
of the local coordinate system, 
(X1,/Y],/2 1); a point on the 
local x @X15 (X9,¥2725) one a 
polnt Im the x->y plane | 
Z5), noncollinear with ene 
first two points and located 
to define desired plus z 


digectron- 
5 OA ee A2 == Ignored 
6 Omit 
di im 7 Card number = 7. 
i 2-4 r3 K Blement number. 0 < Ky )gece 
i 5-7 1 ee 2A Same as Card 1, columns 50/7. 


(may be omitted) 


* s 
femoers given in this column are to be taken literally =e ine 
meer must substitute appropriate numbers for letters. 


me eee 





GABLE so (Conese 


LINE ELEMENT DATA 


Card Columns Format ITtem* Interpretation 


iy Oneal 27 0 SA First and second substitute 
ie Grid point numbers Tet enc 
grid point numbers are nega- 
tive, solution displacements 
are produced in the local 
coordinate system at the 
Grid Powe: 


22-28 i a io fgnored 


PSIG S70 = Coordinates (inches) of the 
substitute grid points corres- 
ponding to the first and 
second elastic grid points: 
all in the local on everaie 
SySten, 


fi au 7 2 A2 & Coordinate identity. If C=L, 
substitute grid point coordi ia 
ates are in the local system. 
If € 4 bL; Coordinates, ame 


overall. 
Omit 
9 al 19 9 Card number = OF 
2-4 re K Element number. QO S Koa oe. 
(may be omitted) 
9 57) dlp, Gor ig’ 2A Same as Card 1, columns 5-7. 
(may be omitted) 
9 8-70 97 20 —-— Coordinates (inches) £662 sene 


fierce (Xy7Y7 721), and second 
(X9,;Vor25) elastic grid points 
of the line element and the 
Gia! poi (X3,Y37,23) defining 
the principal plane of the 
cross-section: all in the over- 
all or local systems; 


9 Le 2 A2 € Coordinate identifications 
C= ib, elastic ged pole 
coordinates are in the local 
system. If C # L, coordinaees 
are in the overall system. 


‘ 
Maumbers given in this column are to be taxen literally |) fae 
Meer mast substitute appropriate numbers for letters, 





Por each material used in Ehe Structuce tvoecards 
of material data are required. These cards give the material 
identification, the temperature associated with the material 
properties, the coefficient of material expansion and 
material stiffness coefficients. The material stiffness 
Seerticients are the coefficients of the stress-strain equa— 


tions in accordance with the following equations: 





The above equations hold true for a monotropic material. If 
Daa = D35 ar Da De = 0 an orthotropic material is des- 


cribed. An isotropic material is described if in addition 


D3 = Des = Deg = Day = Day = Dj and Do4 = (V/s v) (D)4)- 


Then Dee = B/(2(1 + v)) and Di4 = E(L = v7 Gls) Cie oe 


@omwe format is as follows: 





*pesn 
UeSceonE UM Si clao2 el PT PemtOT. Spares ofpaed Telz2oq7eul FO SATeO [Te 4eut SaPoOTPUT OF peSsM 
me =~ ee eee. SO 


00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 COOQ0Q0000 x GuaVvO 


0 9450 °7 0 9450 °7 O+8°9T 0 9+0°S 9+35°9 Z Guwo 
O+0° 0 0 9+5°9T 9+9°8 94+5°9T P-SXSCT° E+0€S° 9LPLOZ T aduwo 
(0° 8d) (0° 8a) KO Sa) (0° 8a) (0° 8H) (O° 6a) (0° 8a) CO= Se) aye oew xc LWWUOT 


(UNUTUNTY otTdoz3OsT) 


DNdNT ATaAVL TYIYALVN WIdNVvs 


‘ut’ bs/sat 199 G 159g SSG 
VVq 1t0q 1¢V¥q «IVGq tsquetotzze09 ssougzjTIs TeTz07eW (0° SAL) ea Z 
: ‘ut’ bs/sct 1€Eq 1CEG 
me eS a guge aoa -SRUSTOTIFSOO SSousgsTAsS Tetso en (0°89) 6-7 ie 
UTHUCY | 
Ssoortbep-yout/seyout :uostuedxes TewArzsy4 jo 4yusTOTFFI9o0O (0° Sa) € ie 
‘soTjyazedorad [TetAsezew usaATtTh ou AOZ sanjerzsdusy uTYyueY (0° 8q) C AE 


“FUP OTT EUOLS Ste 20d 

-“WNU UOTIARPOTITIUSPT AseWZoeAeYyo-XTS SYA Jo sAzszZoOeAeyYoO 

OM} HutpestT euz ATuUO *eTqtssod se AejOsSuT poesn 3q 

STSqWUNU UOTIETOOSSY WNUTUNTY pue gwS pazepueqs ey 

paeoUsSUMIODSeA ST 4I “azequMU UOTIeOTJTJUepPT senbtun e 
SACU FSNU TeTAs_Zew YOeR suOTIeoOTJTIUSpT Tetszozeu ouL Visage = Kc} T T 








EU T Plots pred 


VWiVd DBNdNI SHIGVL TWIYaivn 


=; 


(iia +4 65 504) 
L aTavi 





The material tables must be followed by one card of zeroes 
Pemsignal the end of the material tables. 

the next input 1s the nonzero boundaryecend@eren. 
applied to displacements or loads. Displacements are given 
in inches for the restrained degrees of freedom in the x, y, 
and z directions and in radians for the restrained degrees 
of freedom about the x, y, and z axis. Loads applied in the 
xX, Y, and z directions are given in pounds and those applied 
meet the x, y, and Z axils are given in inch-pounds. ‘The 
loads are those associated with the unrestrained degrees of 
Freedom. Card format is the same as that for the weight 
Mati 1X . 

The last section of input is the title information 
@aecs. Tnere are a total of cleven cords defining the dis- 
placements, generalized stiffness and generalized weight 


Memerreces. This is written as follows: 


mec 1,2,3 = DISPLACEMENT, GENERALIZED STIFFNESS AND GEN- 
BRAL WEIGHT MATRICES 


ROW AND COLUMN CORRESPOND TO MODE, LOAD, OR REACTION 
DISPLACEMENT OR EXCEPT DDC L ROW LEADING DIGITS CORRESPOND 
merGRID POINT AND FINAL DIGIT TO THE DEGREE OF FREEDOM At 
fit GRID POINT 

REACTION NOT DISPLACEMENT AT RESTRAINT 

INTERNAL FORCE FOR ELEMENT -X- (RDC X) 


Row LEADING DIGIUS CORRESPONS TO GRID POINT, FINAL Buin 


1 AXIAL FORCE 4 TORQUE 
8 SHEAR ALONG X2 5 MOMENT ABOUT X2 
7 SHEAR ALONG X3 6 MOMENT ABOUT X3 


COL CORRESPOND TO MODE, LOAD, OR REACT. DISP. OR ACC. 
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The complete computer program is availabie in the 


Naval Architecture and Marine Engineering Library. 





MODELING OF THE STRUCTURE 


The structure analysis of a systemoas Comolijcatced 
as an icebreaker bow requires considerable preparation in 
order to use the SAMIS program. The first requirement is to 
develop a three-dimensional picture of the structure and 
Serent Lt in relation to the common coordinate system. In 
order to clearly indicate the necessary details and label 
the elements and grid points, the scale must be quite large. 
In Appendix A, the three-~dimensional drawing of the ice- 
breaker WESTWIND is ina scale of 1/2" = 1'. Once the 
three-dimensional table has been drawn, a material table 
fierestebe set up to give the necessary material parameters 
required for the input data. Appendix B gives the material 
table for the icebreaker WESTWIND. 

Since a ship's structure is symmetrical about 
the centerline, only one half of the ship need be drawn and 
emiaivyzed. The boundary conditions for the nodes on the 
centerline can be completely described by displacements. 
Unless ena whipoing Of the ship occurs, - all, nodes can 
be fixed in the athwartship direction. Movement in the fore 
and aft direction will be sufficiently modeled by the facets 
edge, and movemant in the vertical direction can be modeled 


as very stiff spring. 


The boundary conditions of the hull are simply the 
average pressure loads for the depth the facet is beneath 
the waterline. At the bow, shell pressure loads for the ice 
contact area can be selected as desired. The total area in 
contact with the ice is that which will give the load res- 
Bouse Getermined by Dr. White (Ref. 5). Shear and bending 
moment diagrams developed from the loading conditions imposed 
can be used to give tne boundary forces where the structure 
1s cut off. These can be applied as concentrated loads in 
the matrix data as was described as Matrix WTR1 in the 
previous section. This will necessitate some hand calcula- 
tions to determine the section modulus and shear area. 
Furthermore, it would be necessary to determine the 
deacceleration at the time of maximum bow load and apply 
that to the mass of the remainder of the hull to geta 
deacceleration force. A check is needed to determine if 
these total forces satisfy equilibrium for the entire 
system being analyzed. A simpler method would be to assign 
large stiffnesses to springs connected to fixed points and 
supply these stiffness as matrix data input as was described 
bieematrix KER] in the previous section. 

It would be desirable to divide the entire section 
being analyzed into small facets and model each stiffener as 
a beam. However, due to the limit on the number of elements 


mie program can handle (999), this type of modeling can only 





be done in those areas where complete and accurate detail 

is desired. The IIT Research Institute Technology Center 
apemicago, Lllinois in LITRI Project J6127 (Rete 19); 
developed what is called an orthotropic plate model for 
stiffened structures so that large facets may be used for the 
remainder of the structure. Thickness is modeled as satis- 


mam the following equation. 


42 2 2 
ae 
where 
aN = bh (the plate area 
Ar = the area of the frame 
2 ae 
_ = Ad a n x + [T + I. 
p ia 


d = distance between the centroids of the plating 
and frame. 


IT. = Moment of inertia of frame alone about its own 
GOncroLdal axis. 


one 
iS 12 Gee 


Db = dietance’ beeween Frames 


h = plate thickness. 


tae modulus of elasticity in the plane of the plate and im 


mnie adirection Enews tLErencres. Bul 1s 





E(A, + A) 
a ee 


9 le 
bt (1l-v =*) 


tiesmodulus of elasticity in the plane of the plate and 


perpendicular to the plate is 


to 
By = iT fy 
and 
= vE 
ay Xx 


The resulting material constants are then 


Diy = y + Be //Ey 
Do, = Fyy + EL/E, 
Do = B, + Be y/By 
D4 = 0 

D235 = Q 

D3, =G 

Dal ~ Bey 

ae ie By 

Dy3 = 0 

Day = BY 

Dee =e 





D = G 
o ox My 
where G = Gi AE 
Utilizing the SAMIS program, version one IITRI analyzed a 
section of the WESTWIND hull from frame seven to thirty-one. 
Results were not compared with stress data accumulated by 
the WESTWIND, but showed a structural weakness at frame 25 


which has been an area of structural failures in the past. 





BUtURE DEVELOPMENTS 


ineecapabelvetres of this program provides ene isi ps 
structure analysts with a very valuable tool. Modeling the 
entire shell of a ship as an unstiffened plate, and then 
designing the stiffeners to be most effective from the 
Seuss OL the first run can lead to Optimal utilization o£ 
material in ship construction not only for tough skinned 
vessels such as an icebreaker, but also for all other types 
of vessels. At last there is a glimmer of hope that we can 
ceparct from tne long standing method of using something tnat 
worked before in a new ship so we can have hove of some 


measure of success. 





HO: 


ee 


We 
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APPENDIX, & 
MATERIAL TABLES 


(All material is medium black steel unless otherwise indicated) 
Saell Plating 

Main deck to second deck. 

20.4# high strength steel plate. 

Second deck to platform deck. 

66.57 high strength steel plate. 

Platform deck to keel. 

Boet+ high strength steel plate. 

Peamaverse Frames 

Seem, Main deck to 30° bow. 

fee 15,.3¢ plate. 

30° bow, frame 7 to keel. 

40" x 20.4% plate. 

ww pow, Lrame 7 to keel. 

20" x 40.8# plate. 

Gant £Lrame 02 to frame 61, main deck to second deck. 
meee”) =X 8 2545 T 

Second deck to keel. 


ieee 4° xX «15.3 T 


Note special framing detail in Appendix C. 





transverse Bulkheads 

(All stiffeners on 24" centers) 
Bulkhead 7 

Pia Giang 

feoot Plate, 

Stiffeners 

oe 4° x Lit T 

Bulkhead 19 

Plating 

Main deck to second deck. 
7.65% plate. 

Beeond deck to platform deck. 
ieee/ 5: Dlate. 

Stiffeners 

oes” 6X 8.254 T 

Bulkhead 31 

Main deck to second deck 
feo plate. 

second deck to platform deck. 
iez+ plate. 

Platform deck to keel. 

feo; Dlate. 

Stiffeners 

Main deck to second deck. 


ex 4° x 1llg T. 


Sep 





Bulkhead 43 

Prating 

Magn aeck to third deck 
fe6ot plate. 

Third deck to keel. 
ier pDlate. 


Seitfeners 


Main deck to second deck. 


Beeex 2.069" x 4.48% T 
Second deck to shell. 
cae 4” CUX:«CLdIF OT 
Bulkhead 61 

Pirating 

Main deck to third deck. 
Weegoy plate. 

Third deck to keel 

ect? plate. 


Stiffeners 


Main deck to second deck. 


pu 2.69" x 4.48% T 
Second deck to shell 
oa 4" x L1F T 
Longitudinal Bulkheads 


meame 31 to Frame oL. 


Tea pe: 





Maan deck to second deck. 

mh. 24 plate 

second deck to shell. 

iz /5+ plate 

stiffeners 

Frame 31 to frame 61 on 16" centers 

coe 4" x L1lft T. 

Main deck 

Plating 

peem tO frame 26 

iz; plate. 

Eee 26 to frame 61. 

feooa;, Dlate. 

Transverse Stiffeners 

Stem to frame 6l 

meee” x 8.25% T. 

meme 4, from centerline to 10' on either side of centerline. 
tae 4" x 16.5¢ T 

Longitudinal stiffeners 

Geame 7 to frame 19, 1' off centerline at frame 7 and 5' off 
centerline at frame 19. 

ee =< 6 1/2" x 25% I 

feenme 31 to frame 61, 9' off centerline at Frame 31 and 19" 


ee CONLerline at Erame 61. 


-~ Fad ~— 





fx) 6 1/2" x 254 1 

Second deck 

Pat aing 

Bow to frame 31 

foot Diate. 

Frame 31 to frame 61 

mic? plate. 

Deel: Stiffeners. All 

Ome 4° X* «Ld TT. 

Transverse beam frame 35 

eee 4" x 25% T. 

Longitudinal beams 

Bow to frame 7 on centerline. 

fee 6 1/2" x 25% I 

imeeecomn Centerline at frame 7 to 6.0 from centerline at 
Lrame 2/7. 

eee > 1/2" x 254 1. 

Frame 31 to Beene £30Om. Centerline’, 
ieee 6 1/2" x 25¢ LT. 

Frame 43 to frame 61, 6.0' from centerline. 
ieee 6 L/2" x 25% I 

Third deck 

Piating (Ali) 

10.2# plate. 

Transverse stiffeners 


Bow to frame ol 
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ee 





fee 6 374" x 15) 1. 

Frame 45 and frame 49 from centerline to 6' on either side 
of centerline. 

ieeex 6 3/4" x 30t T 

Longitudinal stiffeners 

ier rom centerline at frame 7 to 9' from centerline at 
Erame 31. 

Frame 31 to frame 43 on centerline. 

Frame 43 to frame 61, 6' from centerline. 

moe x 7° x 364 I 

Mee7or Drackets, 21" x 10 1/2", used to tie into frame at end 
of each deck stringer. 

Platform Deck 

2 levieskale 

Frame 7 to frame 31 and frame 43 to frame Ol. 

iOmez+ plate. 

Mieansverse stiffeners (All) 

fee 6 3/4" x 30F T. 

Frame 44 and frame 48, from centerline to 6' on either side 
Srmecenterline. 

Meee 6 3/4" x 30F T. 

Longitudinal stiffeners 

Frame 43 to frame 56 on centerline. 


ee x f° X. 30%. TL. 





Deck 


Menglevdinal seh) 


6" x 4" x 16% I 


Equal 






APPENDIX C 


FRAMES 
Shag taal 0. OM EE 





i 2nd Deck : 
| 
l 


Sra beck 


Platform 
Deck 


APPENDIX C 
(Coirmewa) 


All stiffeners spaced 
equa distancenapare 
on 
Shell Frames 
S060 ene a 


oe Gl sc me iat aa 
(ALL) 
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